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In this paper we argue that the Landau singularity of the coupling constant 
in QED reflects the instability of the vacuum state, in other words the pertur- 

VO . 

bative ground state of QED. We do not have a rigorous mathematical proof. 

We rather provide strong arguments derived from the theory of stability of 

matter developed by Lieb et al on the one hand and the theory of weakly 

interacting fermions developed by Feldman et al on the other hand. 
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Landau singularity of the coupling constant in Quantum Electrodynamics has been 
known for more than four decades [1]. The original derivation of the singularity was based 
on the summation of one loop diagrams of vacuum polarization tensor for photons in the 
perturbation theory. In the early days the validity of such an expasion was looked upon 
sceptically by many people including the authours themselves. However, after the advent of 
1/N expansion technique, it was soon realised that this singularity appears at the leading 
order in 1/Nf expansion where Nf is the number flavour or number of specieses of elec- 
trons. This implies that in the infinite flavour limit this singularity is exact. There have been 
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attempts to interpret this singularity in many ways. Landau and Pomeranchuk tried to ar- 
gue that this singularity reflects the fact that at short distances strong vacuum polarization 
effects screen the electric charge completely. Others have called it Landau ghost reflecting 
the internal inconsistency of Quantum Electrodynamics [2]. In this paper, we shall argue 
that Landau singularity reflects the instability of vacuum, in other words the perturbative 
ground state of QED.We do not have a rigorous mathematical proof. We rather rely on 
'analogy and precedence', the only other method available to a theoretical physicist in the 
absence of a mathematical proof.For our 'analogy and precedence' we refer to the theory 
of stability of matter developed by Lieb et al on the one hand and the theory of weakly 
interacting fermions developed by Feldman et al on the other hand. 

In quantum theory a many-body system with ground state energy E is called thermody- 
namically stable or simply stable if E /N is bounded bellow when the number of particles, 
N — > oo. In this context, the question of stability of matter consisting of negatively charged 
electrons and positively charged nuclei is very important.lt was in 1967 when Dyson and 
Lenard proved that, in the framework of nonrelativitic quantum mechanics, matter con- 
sisting of N electrons and K static point nuclei is stable [3]. Subsequently, Lieb and his 
collaborators have made a very detailed investigation of the stability of matter in nonrel- 
ativistic as well as relativistic case ( [4,5] and references there in ). These studies seem to 
suggest that at high enough energies Quantum Eletrodynamics may not be a well behaved 
theory. Here we argue that there exists one approximation scheme for QED in which it is 
really so - the perturbative ground state of QED is unstable. This is the large flavour limit 
of QED, where number of flavours Nf tends to infinity, while e 2 Nf is held fixed , "e" is the 
QED coupling constant. Thermodynamic stability of a system of particles interacting via 
electromagnetic interaction is associated with the control of the short distance behaviour of 
the interaction. In the nonrelativistic limit, zero point kinetic energy of the fermions controls 
this short distance behavoiur- however, in the relativistic limit, there is a need for certain 
bounds on the value of the fine structure constant [4-6]. In the large flavour limit of QED, 
Landau singularity becomes exact and it is well known that this is also associated with the 

2 



short distance behaviour of the electromagnetic interaction [1] . We refer to the analogy of 
this model with the theory of weakly interacting electron system analysed by Feldman et al 
using renormalisation group technique. We find an exact analogy between the Landau singu- 
larity in QED and the singularity of effective coupling constant in some angular mometum 
channel in the condensed matter system. In the condensed matter system, we know that the 
singularity reflects the instability of Landau-Fermi liquid ground (generally called the BCS 
instability). This analogy clearly suggests that the Landau singularity reflects the instability 
of ground state in QED. We shall argue in this paper using the example of large flavour 
QED, that the singularity of effective running coupling constant in the second quantized 
quantum theory is the equivalent of thermodynamic instability in many-body systems of 
particles in the first quantized quantum theory. 

It should be mentioned here that the stability of ground state and the possible existence 
an ultra-violet fixed point has been studied extensively in the lattice formulation of QED in 
the wide range of values of the fermion flavour Nf by Kogut et al [8,9]. On the other hand 
in massless QED in the continuum, Miransky argues that there exists a chiral symmetry 
breaking phase [10] .However, in these studies Landau singularity plays no role. There has 
been some recent studies of Landau singularity using lattice formulation of QED by Gockeler 
et al [11]. These studies seem to suggest that chiral symmetry breaking allows QED to escape 
Landau singularity.But then chiral symmetry breaking, as their study shows, seems to be 
intimately connected with trivility of QED. Our approach is motivated by the problems in 
the theory of stability of matter as the following section would elaborate. We are interested in 
finding the meaning of Landau singularity rather than finding a way to escape it. (However, 
we do believe that it is important to investigate the consequences of escaping the singularity 
as the study in [11] seem to suggest) In our approach Landau singularity is not treated as 
a ghost as in many other publications but as a signature of instability of the underlying 
ground state. 
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I. THE LARGE FLAVOUR QED 



The stability of many-body system of fermions corresponding to large flavour QED can 
be considered in two different ways. One way is to assume that at lower energies only a 
few flavours are excited. In that case it becomes similar to the single flavour case. In this 
way we avoid confusion that may arise due to the double limit, number of particles iV and 
number flavour Nf, both tending to infinity. At high energies when we need to consider 
QED, all specieses of fermions are excited and to avoid additional complication , they can 
be assumed to be of the same mass. However, we believe this is not the correct way- here we 
miss the role played by flavour in a many-body theory. It is possible to consider the quantum 
mechanical many-body problem with large flavours by carefully taking care of the double 
limit mentioned above. Let us consider a region of linear dimension R where there are N 
number of fermions with Nf flavours. Assume N >> Nf. Taking into account the the Pauli 
exclusion principle for the case involoving multiflavour fermions in 3-dimensional space, the 
kinetic energy K can be written [12] as K |^ ^2/3 • The potential energy U due to 

rn R N f 
II 

total energy, E(R), is given by 



Coulomb interaction is given by [4], U — —a a e2N ^ ,i , a > is a numerical constant. The 



h 2 jV 5//3 e 2 A 4 / 3 

E(R) = K + Un- — jj.-av—— (1) 

The ground state is obtained by minimizing the total energy, E(R), with respect to R. 

E - -^{e 2 Nf) 2 N; e = E /N « - ™£{#Nf? 

h 2 N 1/3 1 _ R h 2 1 

When the number of flavour Nf increases, the size of the small box tq decreases and when 
it approaches Compton wave length we can no longer use the nonrelativistic quantum 
mechanics. Therefore, the above estimates are correct only when 



In the relativistic quantum theory of many-body systems, it is almost customary now to take 
c\p\ [4] , where c is the velocity of light, as the kinetic energy of individual particles. Therefore, 
the kinetic energy K for the system of N particles with Nf flavour in a region with linear 
dimension R is given by, K = ^1/3 • Therefore, in the relativistic case, 

„,r,S r, rr ^cN^ G 2 N^ 

E{R) = K + U « "°o— r~ (4) 

In this case stability of ground state requires E(R) > or = O.This leads to the condition, 

N f < or = —77 ^ (5) 

which in the unit c = % = 1 , becomes Nf < or = r^sp • At this ponit a few comments 
are necessary. We considered a theory with large N and large Nf, and it is implicit in our 
discussion that we are considering, N — > 00, but we have said nothing about the flavour 
Nf.We could have also taken Nf — > 00 but slower than N. We could also assume that 
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in the limit Nf — > 00 and e — > , e 2 Nf = constant, and satisfies the equation 
Eq.[5] and conclude that in this limit matter is stable. Mathematically it is correct. But 
this is only a formal result. e does not tend to zero and Nf , never tends to infinity. The 
mathematical limit, has the physical implication that when e is small and Nf is large, the 
relevant parameter that sets the stability criteria is e 2 iV^ 3 .It is the value of e 2 iV^ 3 that 
decides the stability of the system. When e 2 N l J z > , in other words, when the number 

a 

of flavours Nf > 7~2j3 , the many-body system is thermodynamically unstable. 

In Quantum Electrodynamics we find that the relevant parameter is not e 2 iV^ 3 but 
e 2 Nf. In large flavour QED with Nf flavours, when we calculate the connected Green's 
functions using Feynman diagram technique, we find that one loop diagram are of order 
e 2 Nf and higher loop diagrams are of lower order in iV/.This becomes more transparent by 
carrying out expansion in 1/Nf. The essential idea is to replace e— >e/^N~f, and consider 
the Lagrangian, 

C = ]T ¥ (<y ^ + m - ej y/N f 7^) V + -F% (6) 
i=i 4 



where ip is the Dirac field, and and are the potential and the field strength respec- 
tively. In the leading (zeroth) order in 1/Nf , the photon propagator receives radiative 
corrections only from the the single- loop Feynman diagrams. Higher loop contributions are 
of order one or more in 1/Nf. The renormalized running (effective) coupling constant at 
energy, k 2 = A 2 , up to the leading order in 1/Nf , is given by 

2 2 e2 
6e//(A ^ = 1 - (e 2 /37r)/nA 2 /m 2 (7) 

where we have used, e 2 = e 2 (k 2 = m 2 ), for the physical coupling constant. One can easily 
infer that the running coupling constant blows up at, A 2 = m 2 exp(3n/e 2 ) . This is known 
as Landau singularity in QED [1]. It is easy to see that this singularity is associated 
with the short distance (large momentum) behaviour of the interaction. To understand the 
meaning of Landau singularity, here we refer to the renormalisation group analysis of weakly 
interacting Fermi system with short range potential at finite density and zero temparature 
by Feldman et al [7,13]. The iterative renormalization group transformations leads to the 
same type of singularity in some suitably defined running coupling constant. 



II. FELDMAN MODEL OF WEAKLY INTERACTING ELECTRON GAS 

We briefly describe here the Feldman Model of weakly interacting electron gas (see 
reference [13] for details). The model of weakly interacting electron gas studied by Feldman 
et al is a condensed matter Fermi system in thermal equilibrium at some temparature T 
(for simplicity, assume T = 0) and chemical potential /i. On microscopic scale(^ 1CT 8 
cm), it can be described approximately in terms of non-relativistic electrons with short 
range two body interactions. The thermodynamic quantities such as conductivity depend 
only on physical properties of the system at mesoscopic length scale 1CT 4 cm), and 
therefore, are determined from processes involving momenta of the order of ^ around the 
Fermi surface, where the parameter, A >> 1, should be thought of as a ratio of meso- 
to- microscopic length scale. This is generally refered to as the scaling limit (large A, low 
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frequencies) of the system. The most important observation of Feldman et. al. is that in 
the scaling limit, systems of non-relativistic (free) electrons in d spatial dimensions behave 
like a system of multi-flavoured relativistic chiral Dirac fermions in 1 + 1 dimensions. The 
number of flavours iV fa const. A d_1 . It is possible then to set up a renormalization group 
improved perturbation theory in j- around the non-interacting electron gas, where in, the 
large number of flavours N, play an important role in actual calculations. 

A. Free Electron gas and the Multiflavour Relativistic Fermions in 1 + 1 Dimensions 

Let us consider a system of non-relativistic free electrons in d spatial dimensions with 
the Euclidean action, 

S (r,^) = E/ d d+1 xrAx)(id - - /i)^(x) (8) 

The Euclidean free fermion Green's function, G Q aa t(x—y), where a and a' are the spin indices, 
x = (t, x) and y = (s,y), t and s are imaginary times, t > s, is given by, 

/ e -iko{t-s)+ik(x-y) 
(dk)- (9) 

Where we have used (dk) = (2^jTd+rj d d+1 k. In the scaling limit, the leading contributions to 
G Q aa , (x — y) come from modes whose momenta are contained in a shell of thickness ^ 
around the Fermi surface Sp- In order to approximate the Green's function, let us introduce 
the new variables uj, p», po such that kpoo G Sf, Po = k and k = (kp + p\\)oo. If k G S^, 
then p\\ « kp, and we can approximate the integrand of Eq.(4), by dropping p| term in 
the denominator. In other words, 

GlAx-y) = S™> J ^^kP l e^-^G c (t-s,u(x-y) (10) 
where da(u) is the uniform measure on unit sphere and 

fdppdpw e -iw-s)+ip^*-in 

G c {t-s,u(x-y)) = - / " : 11 

J 27r 2n zp — Vpp\\ 
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is the Green's function of chiral Dirac fermion in 1 + 1 dimension, vf = kp/m is the 
Fermi velocity. The cJ-integration in Eq.(5) can be further approximated by replacing it 
with summation over discrete directions ujj by dividing the shell Sp ^ into N small boxes 
Bq v j = 1, .., N of roughly cubical shape. The box, B^, is centered at ujj £ Sf and has an 
approximate side length ^.The number of boxes, N = Jl d _iA d_1 , where is the surface 
volume of unit sphere in d spatial dimensions. The Green's function is, now, given by 

r dr> n dm r> , P -ipo(t-s)+ip(x-y) 

GUx -y) = -<W E / f^^ — (12) 

rr* J 2n 2n 2n ip — v F p\\ 

where p = p\\£u + p± is a vector in B$ i — kpujj and p e 1Z. Thus in the scaling limit, the 
behaviour of a rf-dimensional non-relastivistic free electron gas is described by iV = fi (i _iA a!_1 
flavours of free chiral Dirac fermions in 1 + 1 dimensional space-time. The propagator G c (t — 
s, u(x — y) depends on the flavour index uj .But the energy of an electron or hole with 
momenta k , depends only on py, where p\\ = k.uj — kp and uj — .It is proportional to 
p\\ , just as for relativistic fermions in 1 + 1 dimensions. 

B. Renormalization Group Flow and the BCS Instability 

Large scale behaviour of the weakly interacting system is described by an effective action. 
To see how the effective action is calculated, let us consider a system with Euclidean action 
of the form, 

s(r^) = s (r,^ + Sj(r^) (13) 

where So(ip*, ip) is the quadratic part given in Eq.(3) and Si(ip*,ip) is the quartic interaction 
term. For a weakly interacting electron gas, 

Sj(r^) = gok F ~ d J2 I d d+1 xd d+1 y : ^MM^Hx - y)6(x - y Q WAv)Mv) ( 14 ) 

<j,a' 

The factor k F ~ d ensures that g is dimensionless. The two body potential v(x) is assumed 
to be smooth and short range, and therefore, its Fourier transform v(k) is also smooth. In 



8 



order to calculate the effective actions, we first split the field variable into slow and fast 
modes, ■?/> = ■?/><+■?/>>, where Supp ^> C R x (R d \ S^) (region >) and Supp 
ijj < C R x S^ (region <) . We then integrate out the fast modes using the functional 
integral for fermions, 

e - s <tf = D^D^e- 3 ^^^^ (15) 

Now using the linked cluster theorem, we obtain 

e -s eff (r<^<) = eX p (-S ,<- < Sj > G o +1 < Sr, Sj > G o -1 < 5 j; S T ; 5 7 > G o -...) (16) 

The abbreviations, < a;b > < a;b; c > etc., denote the connected correlators. The subscripts 
G> indicate that the expectations < (.) > G o are calculated using infrared cut-off free propa- 
gators in accordance with the functional measure, dPtyt,,^) = (1 /E.)DipZ > D'ip > e~ So ^>'^ > \ 
The connected correlators can be evaluated using the Feynman diagram technique. There- 
fore, the effective action can be calculated once we know the amplitudes of connected Feyn- 
man diagrams. It is clear that the S e ff contains far more interactions than the original 
quartic interaction Si. However, for weakly interacting systems the original coupling remains 
dominant. To carry out the iterative renormalization group scheme of Feldman et al, we 
choose some large scale Ao << ^= and calculate the effective action, S e ff perturbatively 
to leading order in ^-.The effective action depends on modes corresponding to wave vectors 
located in the shell, S^p \ of width ^ around the Fermi surface. Now, we divide the shell 
Sp^ into iV = const. A d_1 cubical boxes, B^, of approximate side length ^.Next, we rescale 
all the momenta so that, instead of belonging to the boxes B$. they are contained in boxes 
B^ j of side length ss k F . These two steps are generally known as decimation of degrees of 
freedom and rescaling. The renormalization group scheme consists of iteration of these two 
steps. 

Assume that the degrees of freedom corresponding to momenta not lying in Sp^ have 
been integrated out. Let tp a (k),k G R x Sp^, denote these modes. The sector fields are 
defined as 
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^{x) = I {dk)e^ kot -^- kF ^ s ^ a (k) (17) 

It is easy to see that ip a (x) = e lkFLJX ^,a(x) . Inserting the Fourier transform of sector 
fields in Eq.(3) and Eq.(9), and carrying out some algebraic manipulations, we obtain 

So = -E[ , , Jdp)^>)(iPo-vpa^ + 0(^))^(p) (18) 
Sj=f4~ d £ vikpfa-X*)) J{dp l )...{dp A ){2n) d+l 5{p l +p 2 -p^-p A ) 

L0i,..,W4;cr,<T' 

^kAP 1 ^k,a'(P2)^ 3 ,a'(P3)^APd + terms of higher order in^ (19) 

Using cluster expansions to integrate out the degrees of freedom corresponding to momenta 
outside the shell Sp^°\ one can show that , at scale ^ , the effective action has the form 
given by Eq.(13) and Eq.(14) except that v(k F (ui — cD 4 )) is replaced by a coupling function 
g(ui, lj 2 , £3, uj a ) ~ v(kp(uji — £u A )) with u)i + u 2 = uj 3 + u A . 

Next, one considers the rescaling of the fields and the action. The fields are rescaled in such 
a that the supports of the Fourier transformed "sector fields" are boxes, = \(Bj — k F u) 
,of roughly cubical shape with sides of length k F , and such that the quadratic part of the 
action remains unchanged to leading order in ^-The first condition implies that p 1 — ► p = p\ 
and x 1 — > £ — f ■ The rescaled sector fields and their Fourier transforms are given by, 

^, CT (0=A"VMAO ; h,AP) = ^ d ^Aj) (20) 

Inserting the scaled Fourier transformed fields into quadratic part of action So, it is easy to 
see that So remains unchanged if the scaling dimension a — |. Now, inserting the rescaled 
fields in the quartic part of the action, we find that the quartic part has scaling dimension 
(1 — d). The quartic terms of higher degree in momenta as well as terms of higher degree in 
fields appearing in the effective action have smaller scaling dimensions. 
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Thus the effective action in terms of scaled sector fields is, 

1 1 



( rf P)^, CT (p)(iPo - v F o;p)^ i(T (p) + -^rj g{O u u 2 ,u 3 ,u 4 ) 
(dpi). . . . (dp 4 )(27r) d+1 5(pi +p 2 - p 3 ~ ^^^(pO^^fe)^^^)^^^) 



+terms of higher order in— (21) 

A 

We see that that the inverse propagator for the sector field is diagonal in uj ,and it depends 
only on p Q and p\\ = up but not on p± = p— (uj.pjuj. 

We are interested in the renormalization group flow equations in the leading order in 
j. Therefore, for carrying out the decimation of degrees of freedom, we will be interested 
only in those diagrams that contribute to the amplitude in the leading order in ^.Before we 
consider these diagrams, let us consider the possible inter sector scattering geometries. How 
many independent g^(ui, u) 2 , uj 3 , 0J4) exists? For d — 3, suppose cu 3 7^ — cJ 4 .On the unit 
sphere, there are = Const. Aq _1 different cD's. But all choises of Q\ and uo 2 with uji+uj 2 = 
uj 3 + C04 lie on a cone containing u 3 and cJ 4 with u 3 + d3 4 as the symmetry axis. Therefore, 
there are 0(Aq~ 2 ) choices. Only when oo 3 = — cJ 4 that there are = Const. choices. 
Couplings involving incoming states with uj 3 7^ — u; 4 will be represented by g^(u)i, uj 2 , uj 3 , cD 4 ). 
Couplings that involve sectors u 3 = — cJ 4 or equivalently Cui = —uj 2 will be denoted by 
g^g CS {uj 1,004). Because of rotational invariance, g^ cs {Q\,(jJi) is a function of only the angle 
between oj\ and cJ 4 . 

The chemical potential receives corrections from connected diagrams with two exter- 
nal electron lines (self energy correction of the electrons). The contribution of order zero 
in j comes only from the tadpole and turtle diagrams. These diagrams contain one in- 
ternal interaction squiggle of order -^rr and there are = Const. Aq _1 choices of the 

A o 

inner particle sector. Therefore, the contribution is of order zero in j-. The amplitude 
corresponding to these tadpole and turtle diagrams turns out to be ^-independent but 
of order O(g^0)).Thus, S/ii = O(g^/X ) , and the renormalized electron propagator is 
Gf (Po,P-v) = - [ipo - v F p.w + A 5/ii] _1 . 
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To find the evolution of the coupling constant g(ui, 002,^3,^4), we have to calculate 
amplitude of diagrams with four external legs. It is found that when Qi + uj 2 = ^3 + ^4 7^ 0, 
the coupling functions do not flow in the leading order in j. But for sector indices uji + uj 2 = 
£3 + 004 = 0, the coupling functions, g^ cs {uji, 0J4), flow. The diagrams that contribute to the 
flow equation are the ladder diagrams with self energy insertion for the internal electron 
lines but with no other two legged subdiagram.The amplitude of such a diagram with n 
interaction squiggles and with zero incoming and outgoing box momenta of the particles is 
given by, 

("^T +1 E {-^) n (3 n gBCs{u,UJ n )g B Cs{tin,tin-l)---gBCs{£l,u) 

W\,...,Wn 

In the equation above, (a;', —a;') and (u;, — u) are sector indices of incoming and outgoing 
electron lines respectively. Other sector indices correspond to the internal electron lines. (3 is 
a strictly positive number coming from the fermion loop integration in the Feynman diagram 
and is given by f3 — J dk±dk\\dk [k% + (v F k\\ — A<5/zi) 2 ] -1 . We find that the renormalized 
value of Qbcs is 

g&M*) = gfcs^^') + E (^) n E (-i) B ^5^,^)&(^»,^i)-^5(^, 

n=l aTi,...,0M 

+o( 9 —) 

The explicit expression for flow equation can be obtained by expanding the coupling func- 
tions g B cs(u, v') = gBcs((u, u') into spherical harmonics, g B cs((u, = £ 9ihi({fi, £'). Up 
to terms of order ^, 

^ +1) = r4^) + °4) (23) 

To obtain the differential equation for the R.G. flow, let us define := g\ Xl \ and consider 
a scale A = e*A . Next, define gi{t):—g\ e X °\ The coefficient (5 = fl(t,t') vanishes in the limit 
t' \ t, therefore, /3(t',t) = (t' - t)-f(t) + 0((t' - tf). Writing difference equation for the 
couplings g\^ +V> and and dividing both sides of the difference equation by (t' — t), we 
finally obtain 
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j t 9i{t) = -l 9l m) 2 + 0{e- t g{tf) (24) 

where 7 = j(t) > . It is independent of / and approximately independent of t , and 
therefore, we set 7 = 7 . The positivity of 7 follows from slow monotone growth of f3(t' — t) 
in t'. Neglecting the error term, the solution can be written as, 

9i(t) = 9l{0 \ , • (25) 

If there is an angular momentum channel, / , with attractive interactions (57(0) < 0) the 
flow diverges at a finite value of the scaling parameter, t = —('Jogi(^))~ 1 ■ This singularity 
reflects the BCS-instability of the ground state. 

III. INSTABILITY OF THE PERTURB ATIVE GROUND STATE IN QED 

We have seen that in the case of large flavour QED as well as the weakly interacting Fermi 
system with attractive potential in any one of the angular momentum channel, the running 
coupling constant blow up at some finite value of their arguments. In the latter case, we know 
that the perturbative ground state is unstable (BCS instability): perturbation theory is built 
around the wrong ground state (the non-interacting Fermi gas). The elementary excitations 
are the electrons and holes, and Fock space is constructed using these single particle states. 
The BCS instability implies that such excitations do not exist- we require a modification of 
the excitation spectrum. But this also means that the new stable ground state is not in the 
eigenstate of particle number operator, N and the global U(l)- gauge symmetry associated 
with the conservation of particle number (excitations) is broken. The BCS ground state of 
superconductivity [13] contains all these features. 

Thus the BCS theory of superconductivity contains two parts : 

1) The Landau-Fermi liquid ground state is unstable 

2) There exists another ground state- superconducting state with Cooper pairs as the 
elementary modes. 
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In this paper, we have been concerned with the first part of the problem. We have 
seen that , up to leading order in 1/N, interactions in the various angular momentum 
channels decouple, and attractive interaction in any angular momentum channel leads to the 
instability of (Landau- Fermi liquid) ground state. This analysis in the Feldman theory was 
carried out using renormalization group technique in the language of " Second Quantization" . 
On the other hand the Lieb theory of stability of matter is discussed in language of "First 
Quantization" . Both the methods address the issue of stability of " Ground State of a Many- 
Body System" . To understand the connection between the two methods, let us consider a 
realistic Fermi system. The most interesting realistic Fermi system happens to be the liquid 
He 3 [14] . The minimum of the interatomic potential in He 3 (as well as in He 4 ) is too close to 
the atoms and is too shallow and, therefore the zero point is high compared to the potential 
energy and the helium atoms do not bind to form molecules. However, at higher density (at 
high pressure) the net potential energy decreases as —N 2 v (v being the two body potential) 
and overcomes the zero point energy which increases as N 5 ^ 3 with the particle number N. 
Thus at higher density, we have the atomic helium liquid which is a Landau-Fermi liquid. 
As the density is increased further potential energy further decreases. At very high density 
the zero point energy is completely overcome, however the repulsive hard core (two body) 
potential comes into play.In He 3 , in the angular mometum channel / = 0, the repulson is 
too strong. However, in I — 1 and I = 2 channel the attractive potential energy over takes 
the repulson [14-16]. Thus the Lieb theory leads to 'unbounded from bellow' energy in 
the two channels (interaction in different angular momentum channel decouple). This makes 
the quantum mechanical Landau-Fermi liquid ground state unstable. On the other hand, 
Feldman theory implies that for the angular momentum channels I = 1 and 1 = 2, the 
effective running coupling constants blows up. 

In the large flavour limit, the effective coupling constant in QED blows up exactly like 
the effective coupling constant in the Feldman theory. The Lieb treatment of the Coulomb 
system also suggests that in the relativistic limit for some value of the number of flavour Nf, 
the ground state become unstable. These similarities clearly suggest that the perturbative 
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vacuum state of QED is unstable. Note that in the Lieb theory for coulomb systems, we 
had nonrelativistic treatment followed by a relativistic system for the stability of the ground 
state. In condensed matter systems, the formation of cooper pairs changes the situation dras- 
tically. The liquid makes a transition to superfluid phase. The cooper pairs are essentionally 
free and therefore, the new superfluid state does not suffer from the type of instability men- 
tioned above. Does the instability of ground state in QED imply the existense of BCS type 
of ground state? The answer seems to be negative. The residue at the poles of the effective 
running coupling constant in QED and in the condensed matter have opposite signs. In the 
condensed matter system it signifies the presence of cooper pair while in the case of QED 
the opposite sign signifies the absence of such pairs [2,17]. 

The remarkable analogy of the evolution equations of the effective coupling constants in 
large flavour QED ( Eq.(2) ) and the condensed matter system above ( Eq.(20) ), lead to 
the suggestion that the perturbative ground state (vacuum state) is unstable in large flavour 
QED too. Landau singularity, just as in case of the condensed matter model, reflects the 
fact that we have constructed the perturbative theory around the "wrong" unstable ground 
state. In our model of QED, the elementary excitations around the perturbative vacuum or 
ground state are the electrons, positrons and photons and, therefore, instability of ground 
state or the vacuum implies that such excitations have become unstable. The new stable 
ground state would require modification of the particle spectrum. There are suggestions in 
the literature that the electric charge charge is totally screened by the new nonperturbative 
ground state [1,8,9,11,18] . The most interesting and exciting possibility is the one suggested 
by Witten [19]. His conjecture is that at very high energy the states carrying electric charge 
become very massive and are suppressed while magnetic monopoles become lighter and are 
excited. 
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